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Abstract: Consequences of a new force mediated by a light scalar particle for neutrino 
oscillation experiments are considered. Such a force could give rise to neutrino masses 
and mixings whose matter dependence is much more significant than the MSW effect. We 
consider the effects of such a new force on the limits derived from oscillation experiments, 
and examine how the constraints on neutrino models are altered. Re-analysis of neutrino 
data, as well as new experiments in which large matter effects are systematically explored, 
is required to directly probe such physics beyond the standard model. 
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1. Introduction and Motivation 

The last ten years have brought great progress with the first direct experimental evidence 
for the physics beyond the standard model, the discovery of neutrino masses and mixing. In 
the standard model, the neutrino is a massless particle with matter interactions mediated 
by W and Z bosons. Wolfenstein [1], Mikheev and Smirnov [2] were the first to show 
that these standard model interactions give rise to matter dependent masses and mixings 
from coherent forward scattering in the medium, and that these effects can be significant 
despite the small cross-sections for such interactions. Subsequently, many authors have 
considered the medium effects of beyond the standard model neutrino interactions, with 
an emphasis on flavor changing neutral interactions (for a recent analysis, see [3, 4] and 
references therein) . These medium effects arise from new interactions whose energy scale is 
above the standard model, and have been shown to be a small effect. The effects of a new 
scalar particle coupling only to neutrinos [5], neutrino mass arising from matter-neutrino 
interactions through such a scalar [6], and the matter effects of a massless string dilaton 
field coupled to neutrinos [7] have also been considered. 

In this paper, we wish to consider matter effects on neutrinos due to interactions with 
a new, very light and weakly coupled scalar particle. In [8], such a new light particle 
was introduced as the connection between neutrinos and the dark energy, to explain the 
apparent coincidence that the dark energy scale (10 -3 eV) 4 is of about the same scale as 
the neutrino masses 5m u ~ O(10~ 2 eV). It was also shown [9] that such new interactions 
with a light scalar particle could give rise to masses and mixings which are environment 
dependent. Here we parameterize such a new interaction by Yukawa couplings, 
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V = \ N <f>NN + Y,K<t>Wi + M<t>), (i-i) 

i 

where <p is the new scalar field, N the nucleon field, and are neutrino fields, written here 
in a basis where interaction with the scalar field is diagonal. A light scalar field coupling 
to neutrinos has been considered before, but without coupling to matter. In a mean field 
approximation, such an interaction would give rise to a neutrino matter potential, which 
is an effective shift in the neutrino mass, 

M?> = ±£i„ N , (1.2) 

m (j> 

where is the nucleon number density, = V^" (</>), and 4> is the vev of 4>, the solution 
of V'((J)) = 0. This minimization equation is equivalent in the mean field approximation to 

Ajvnjv + K n »i + V o(4>) = 0- (1-3) 

i 

Now if Vq is a sufficiently flat potential, which is the case if (f) is a light particle, then 
<fi will vary significantly with (and n v ., which is not of interest for the environments 
considered here). This implies that there are two sources for density dependence in : 
not only the factor of in MfJ^ , but also the dependence m^(<^(njv))- 

M^* in turn modifies the effective neutrino Hamiltonian H e ff, which in a two neutrino 
formalism reads, 



^ e// " 2E U { M^ 2 (m 2 + M 2 ^)^ J ^ + ^ J ' (L4) 

where U is the mixing matrix, \^2Gfn e is the usual MSW term, m; are the vacuum neutrino 
masses, and M- is some linear combination of those given in eqn. (1.2). The result of 
these new terms in H e ff is that the effective mass splittings and mixings (determined 
from the eigenvalues and eigenvectors of the matrix), and their density dependences, are 
modified. The modification could be quite significant, as the new terms could yield effects 
potentially much larger than the standard MSW term. Consistency with gravitational 
force experiments requires the coupling in the new term A at < 10~ 2 rriN/Mpi [10]. This 
small coupling, however, can be compensated by tin in typical earth densities 3g/cm 3 , and 
would give rise to a potential of size 

where p° N = 3g/cm 3 , so that MjJ^ could potentially be much larger than the MSW term, 
of O(10 -23 eV) for typical earth densities. 

In this paper we will take the point of view that these non-standard matter effects do 
in fact dominate over the MSW effect. We will then consider how our interpretation of 
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neutrino oscillation experiments shifts within the framework of these new matter effects. I 
wish to stress that, although such a new light particle was originally motivated in connection 
with the dark energy, non-standard interactions with any new light scalar particle, due to 
the long range nature of the force enhanced by n^, could potentially give rise to new 
significant matter effects. This is in contrast to the effects considered from interactions 
with new heavy particles, which are small due to the short range nature of the force. 

The outline of the paper is as follows. In the first section we make model independent 
considerations on how the interpretation of neutrino oscillation data is changed when large 
non-standard interactions are considered. In the second section we apply these new con- 
siderations to the standard three neutrino model, and discuss the implications for future 
experiments. Lastly, we specifically consider the range of possible models consistent with 
results from LSND. 

2. New Limits from Medium Dependence 

For simplicity in the analysis, and because we wish to make our considerations independent 
of a model for how M- depends on density, we consider here only two media: air versus 
rock, or more generally, any High Density Material (HDM hereafter), which includes any 
combination of rock, steel, and concrete. We assume that the effects of rock or any HDM 
on the oscillation parameters are similar, and that they are significantly different from 
the air parameters. Additionally, because the analysis here is independent of a model 
for how M ? 6 depends on density, no limits from solar neutrinos are derived. To extract 
the vacuum neutrino masses and mixings from any solar neutrino data requires a model 
for how H e ff and hence m^(</>) depends on density; whether the masses and mixing are 
significantly different from the terrestrial values is highly model dependent, and hence will 
not be considered here. Lastly, we assume that at the interface of two different media the 
sudden approximation is valid, so that the evolution of the neutrino phase is continuous. 
This is valid so long as the Compton wavelength of the neutrino mass in the medium in 
short compared to the distance traveled in the medium. 

With these assumptions, we summarize in table 1 the experimental signals we have to 
date, their HDM and air pathlengths, and limits derived. There are several modifications 
to the limits derived when considering large, non-standard in-medium effects. First, as 
mentioned above, in order to derive limits from the solar neutrino data, we need a model 
to determine how Vq'(0), and hence M- , depends on density. Without such a model, we 
are unable to extract 5m 2 and sin 2 29 from the solar neutrino data, as the non-standard 
dependences of <5m 2 and sin 2 29 on solar density are unknown. 

Second, the standard analysis of the Super-K data utilizes an up-down ratio of v e and 
events [11], one pathlength of which is mostly in air, and the other pathlength mostly in 
rock, which could potentially have very different mixings and splittings. To find the range 
of allowed parameters in HDM, one can instead make use of the ratio of upward stopping 
to upward through-going muon sample [12], with limits as given in the final column of the 
table. Notice that these limits are considerably less restrictive than those of the whole set. 
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Third, for negative searches partially in rock or air, Am^ in is scaled by the fraction 
of the path in HDM. One can see this by noting that for small Am 2 L/AE, the amplitudes 
are simply additive: 

A(v x -> Vy) = {v x \U n ...U 2 Ui\v y ) = A n (y x -> v y ) + ... + A 2 (y x -> v y ) + A 1 (v x -> v y ), (2.1) 

where U n is the evolution matrix in medium n, and A n (v x — ► z^ y ) is the oscillation amplitude 
in medium n. To determine the bounds in the medium of interest, we set all other oscillation 
amplitudes to zero except in the medium of interest, and we find that in medium n 

P n (y x - vy) ^ (Aml m L n /2E) 2 , (2.2) 

so that Am 2 min is simply scaled by the ratio of the pathlength in medium n, L n , to the 
total neutrino pathlength. To determine sin 2 20 m i n from these null searches, we again take 
mixings in all other media except the one of interest to zero. Hence, at large Am 2 , 

P(v x -> v y ) = - sin 2 26 min , (2.3) 

so that sin 2 29 m i n is unchanged. This applies, for example, to the CHOOZ reactor ex- 
periment, where the upper bound on Am 2 in HDM is scaled by a factor of 5-10, as the 
experiment is 80 — 90% in air, but at large Am 2 , sin 2 29 m i n = 0.1 is unchanged. Therefore, 
the limit on v e disappearance in HDM for small Am 2 is derived, not from CHOOZ, but 
from the Palo Verde reactor experiment (« 95% in HDM). 

Finally, we note that we have no signals, and only constraints from null searches, for 
the neutrino parameters in air. Therefore, in what remains in this paper we will consider 
only the range of possible neutrino models admitted by the data in HDM. We will find that 
the range of possible models and allowed parameters is quite different from those admitted 
in the standard picture with standard interactions. 

3. Three Neutrino Models 

Within a three neutrino framework, the experimental signals given in the table are tradi- 
tionally understood as follows. The parameter space is determined by two mass splittings 
and three mixing angles. One mixing and mass splitting is determined from the solar 
neutrino data, which is interpreted as the MSW conversion of v e — > v^, v T in the sun, 
with large mixing and mass splitting O(10~ 5 eV 2 ). The other mass splitting and a second 
mixing angle is determined from Super-K atmospheric data, which is interpreted as the 
vacuum oscillation of v^ — > v T , with large mixing and mass splitting O(10~ 3 eV 2 ). The 
limit on the third mixing angle comes from the CHOOZ experiment sin 2 26\^ < 0.1. The 
combination of these limits is illustrated in figure la. 

The constraints, however, are significantly modified with the new interactions consid- 
ered here in HDM. This is illustrated in figure lb. The limits on v e disappearance in HDM 
are derived from KamLAND, with the approximate allowed range 10" 5 eV 2 < Am 2 < 
10~ 3 eV 2 , which is much less restrictive the limits from SNO 6.5 x 10~ 5 eV 2 < Am 2 < 
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Table 1: Range of allowed mass splittings for standard interactions (SI) versus large, medium 
dependent interactions. Am^ in max defines the range of values allowed by the signal at 90% C.L. 
Am^„ is the upper bound on Am 2 at sin 2 29 = 1 at 90% C.L. 



8.3 x 10 -5 eV 2 . For disappearance, the Super-K upward muon data constrains the pa- 
rameters in HDM, which is less restrictive in both Am 2 and sin 2 29 constraints than the 
full set utilizing an up-down ratio of events. But although we cite the upward muon data 
set as the best limits in HDM on disappearance, bear in mind that even these limits may 
not be a rigorous measurement of the HDM parameters, as the analysis contains contami- 
nation from a horizon bin whose pathlength is mostly air. For a truly rigorous constraint, 
a reanalysis of the data would have to be conducted where the horizon bin is removed. 

These modifications in constraints have several implications for the interpretation of 
the experiments in a three neutrino model. In a conventional analysis CHOOZ limits 
sin 2 26>i3 < 0.1. With the new interactions considered here, however, the most obvious 
limit on sin 2 2#i3 in HDM is derived from Palo Verde. For the allowed region of Super-K 
parameter space, Am 2 > 1.5 x 10 -3 , it requires sin 2 2#i3 < 0.5, a rather unrestrictive limit. 
In addition, the Super-K experiment has not carried out an analysis on v e disappearance 
mixing angles with only rock pathlengths to restrict sin 2 2#i3 . The limit that we do have 
on atmospheric v e disappearance was derived in the full set utilizing an up-down ratio of 
events, with both air and rock pathlengths, in order to remove uncertainties as large as 20% 
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Figure 1: Allowed regions for — > i/ x and v e — > z/ x for a) standard interactions and b) non- 
standard medium dependent interactions in HDM. The curves are parameterizations of the pub- 
lished allowed regions. 

in the atmospheric neutrino flux. For z/ e , this-up down asymmetry A = (U — D)/(U + D), 
where U is the number of upward-going events and D is the number of downward-going 
events, is measured to be quite small A Ve = —0.009 ± 0.042 ± 0.005 [33], indicating that to 
high certainty the observed and produced u e fluxes are the same. 

To derive, however, the allowed mixings in air and rock from A Ue and A Vjl would require 
a new fit with six free parameters, 

Am 2 HDM , sin 2 2#f 3 DM , sm 2 29g DM , Am 2 air , sin 2 2^, sin 2 20$, (3.1) 

instead of two. As will be shown here, even an uncertainty of 5% in the ratio of observed to 
produced v e fluxes for maximally oscillated neutrinos could potentially allow sin 2 2#i3 to be 
quite large. To quantify how large #13 can be consistent with Super-K's lack of observation 
of v e oscillation, we write the ratio of the observed v e flux to the v e flux produced in the 
atmosphere: 

We consider the limits from a maximally oscillated neutrino, which gives the tightest con- 
straints on sin 2 2#i3 . Most of such neutrinos have energy ~ 1 GeV and are coming directly 
from the opposite side of the earth (cos 6 ~ — 1). For these neutrinos, the ratio of to v e 
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flux is [34] 

W - 2 - 5 > (3-3) 

and is quite well known (uncertainty ~ 5%). Substituting this ratio into the relevant 
oscillation probabilities, we find 

5p = 1 - sin 2 20i3(2.5 sin 2 6 23 - 1) sin 2 x, (3.4) 

where we are using the usual conventions #23 for the fi— r mixing and #13 for e-disappearance, 
and x = Am^ tm L / AE . To get a rough estimate of the limits on #13 in HDM from Super-K, 
we note again that v e are not observed to oscillate within some fractional uncertainty e in 
flux, so we can approximate 

l-e<|^<l + e. (3.5) 

For neutrinos which have maximally oscillated, sin 2 x = 1, we find 

-e < sin 2 26»i 3 (2.5 sin 2 6> 23 - 1) < e. (3.6) 
Likewise, we can write for disappearance 

= _ ^) + _Lp(^ e _ (3.7) 

so that 

-j^ = 1 - (2.4 sin 2 (9 23 cos 2 6>i 3 - 4 sin 4 (9 23 cos 4 9 13 + 1.6 sin 2 6> 23 cos 4 6 13 ) sin 2 z. (3.8) 

The Super-K upward muon analysis requires sin 2 1Q a tm > 0.8. Then we have the second 
constraint 

0.8 < 2.4 sin 2 (9 23 cos 2 13 - 4 sin 4 (9 23 cos 4 13 + 1.6 sin 2 <9 23 cos 4 6 13 < 1.0. (3.9) 



Shown in figure 2 is the allowed regions given by eqns. |3.6| , |3.9| for e = 0.05,0.1 and 
0.2. The overlap of the constraints gives the combined allowed region. We can see that 
while sin 2 2# 2 3 is constrained to be quite large, there is much less constraint on sin 2 2#i3 
consistent with the data. In particular, for sin 2 2823 = 1> srn2 2#i3 could be as large as 0.2 
for e = 0.05, or as large as 0.8 for e = 0.2. Such solutions can occur when v e depletion is 
canceled by replenishment through oscillation to within an amount e. We conclude that 
sin 2 20i3 could be much larger than concluded from the standard analysis. 

A large sin 2 2#i3 would also indicate that KamLAND could observe very significant 
v e oscillation to the Super-K mass splitting. This possibility is excluded in the standard 
three neutrino model through the limits of the CHOOZ experiment, an experiment mostly 
in air, sin 2 2#i3 < 0.1. As we have shown here, however, the new best limits from Super-K, 
KamLAND, and Palo Verde cannot exclude this possibility in HDM, and in this case a 
three neutrino analysis of the KamLAND data would be necessary. 
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Figure 2: The curves marked e = 0.05,0.1,0.2 are the constraints from v e disappearance from 
Super-K given by eqn. (3.6) for 5%, 10% and 20% uncertainty in the neutrino flux. The upper 
curve is the constraint for sin 2 29 atm from eqn. (3.9). The shaded regions are the overlap of these 
constraints which gives the allowed values of sin 2 2#23 and sin 2 29\3 from Super-K. Dark shaded is 
for e = 0.05, medium shaded for e = 0.1, and light shaded for e = 0.2. 

To this point we have been assuming that Super-K and KamLAND are observing 
separate mass splittings, as the Super-K best fit point lies around 3 x 10 -3 eV 2 and the 
KamLAND best fit around 5 x 1(T 5 eV 2 . But we also note that the new limits from 
Super-K upward muon set on v„ disappearance and the KamLAND data nearly overlap 
in Am 2 . One could ask then whether the two signals could be accommodated with a 
single mass splitting in HDM. Now KamLAND reports a signal with sin 2 26 ex > 0.4. Such 
a large mixing for v e disappearance would have to consistent with the Super-K lack of 
observation of v e disappearance. As we have shown, sin 2 2#i3 = sin 2 28 ex > 0.4 is allowed 
for e > 0.1, which may very well be consistent with the Super-K result. We conclude that we 
cannot rule out the possibility that KamLAND and Super-K are observing the same mass 
splitting. But again we note that although it appears that the fits of both experiments can 
be accommodated with a single splitting, it is only accommodated through less preferred 
values of the mass splitting, as the Super-K upward muon best fit lies around 3 x 10 -3 eV 2 
and the KamLAND best fit around 5 x 10~ 5 eV 2 . 

To narrow these allowed regions, a reanalysis of the full Super-K data set would need 
to be carried out allowing six free parameters on and u e disappearance, as listed in 
eqn. (3.1). Once this analysis has been carried out, one can determine whether the com- 
bination of these parameters with the Palo Verde limits constrain #13 such that a three 
neutrino analysis of KamLAND is necessary. In addition, in the near future we will have ad- 
ditional signals when the MINOS and OPERA experiments measure Am 2 HDM , sin 2 26^ M 
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through direct detection of the appearance of Vn — ► v T . MINOS will also be able to probe 
sin 2 26^ DM through — > z/ e down to mixing angles sin 2 26± 3 DM ~ 0.1. We have shown 
here that, with these new matter effects, the value of sin 2 26^ 3 DM could potentially be 
much larger than is now thought. Should MINOS find a value of sin 2 29\s hi HDM which is 
larger than allowed by CHOOZ experiment, we have a smoking gun for new matter effects 
in three neutrino models. 

4. Models Consistent with LSND 

Next we consider how these new matter effects improve the fit of 3+1, 2+2 and three 
neutrino models to the data with LSND. Shown in fig. 3a is the LSND allowed region for 
Vn v e oscillations, along with the best constraints for v e — > v x and v„ — > u x oscillations 
from Bugey and CDHS. The modified constraints in HDM are shown in fig. 3b. The 
Palo Verde experiment, not Bugey, constrains sin 2 29 e in HDM. The limits from CDHS 
on disappearance in HDM are unknown. 1 In addition, the KARMEN experiment, 
which is ~ 50% air no longer constrains LSND. We consider the implications of these new 
limits in HDM on 3+1 and 2+2 models with LSND, which are disfavored with standard 
interactions. We also consider how a three neutrino model with LSND becomes workable 
with non-standard interactions in HDM. 

4.1 3+1 models 

These medium effects improve the fit of 3+1 models to the KamLAND, Super-K, Bugey, 
Palo Verde, and LSND results. 

In typical 3+1 models, we have three splittings, O(10~ 5 eV 2 ), O(10~ 3 eV 2 ), and the 
larger LSND splitting 0(0.1 — 1 eV 2 ). With only standard interactions, 3+1 models are 
disfavored [35], and a 3+2 scenario [36] improves the fit somewhat, but only by contriving 
to set the masses of the sterile neutrino in the region where the constraints from CDHS 
and Bugey are weakest. 

The reason that these models are disfavored is that the combination of the Bugey and 
CDHS limits on v e and result in a \i — e mixing angle which is too small at 90% C.L. 
to accommodate the LSND result. To quantify this observation, we write the v e and 
disappearance angles measured by Bugey and CDHS in terms of the fraction of v e or in 
the LSND mass splitting, U e4 and U^: 

sin 2 20 e = 4(1- | U e4 | 2 ) | U e4 | 2 (4.1) 

sin 2 2^ = 4(1- | UyA | 2 ) | | 2 • (4.2) 
Likewise, the LSND jjl — e mixing angle is written in terms of these same quantities 

sin 2 20^ = 4 | U e4 | 2 | | 2 = (1 - cos 20„)(1 - cos 20 e ). (4.3) 

x My thanks to Joe Rothberg and his collaborators at CERN for looking into beam path information. 
The lapse of time between the experiment and the present, however, has not allowed them to track this 
information down. 



- 9 - 




Figure 3: LSND allowed region along with constraints on v e and oscillations from null searches 
a) Bugey and CDHS with standard interactions and b) Palo Verde and CDHS with non-standard in- 
medium interactions. As the fraction of the CDHS bcamlinc in HDM is unknown, the limits for non- 
standard interactions are at most as restrictive as those shown. The curves are parameterizations 
of the published allowed regions; in addition irregularities in the Bugey curve for large Am 2 have 
been suppressed, as they are irrelevant for the analysis conducted here. Dashed lines indicate limits 
from null searches, and solid lines allowed regions from positive signals. 

When the limits from Bugey and CDHS on sin 2 29 e and sin 2 29^, as shown in fig. 3b, are 
substituted into these expressions, we find that the resulting allowed mixing sin 2 29^ is too 
small at 90% C.L. to accommodate the LSND result. There is only a small region of LSND 
parameter space, 0.2 < Am 2 < 0.3, which, although excluded at 90% C.L. is allowed at 
99% C.L. Hence these fits are only marginally allowed within a small region of parameter 
space. 

We immediately see how the situation changes with these new matter effects. Bugey, 
an experiment mostly in air, can no longer constrain LSND. Instead, the constraint on 
v e disappearance in HDM comes from Palo Verde sin 2 29 e < 0.17, which corresponds to 
the region of LSND parameter space Am 2 > O.leV 2 . These limits are considerably less 
restrictive than sin 2 29 e < 0.04 from Bugey, and the result is that the combined limits of 
Palo Verde and CDHS do not disfavor the LSND result. 

The Palo Verde limit implies that if MiniBooNE sees e — fx oscillations, it should see 
them with Am 2 > O.leV 2 . Furthermore, if oscillations are detected at MiniBooNE in the 
region O.leV 2 < Am 2 < 0.25eV 2 (excluded by Bu gey in a conventional analysis), we have 
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strong indication for non-standard matter interactions for neutrinos. 

An additional prediction of this model is that a reactor experiment, with similar sensi- 
tivity as Bugey, but with neutrino pathlengths instead in HDM, would detect v e disappear- 
ance, provided the CDHS limits are rigorous in rock. The reason is that the combination 
of CDHS and Bugey limits with standard interactions disfavor at 90% C.L. mixings large 
enough to accommodate the LSND result. We would expect a signal in the approximate 
range 0.04 < sin 2 29 < 0.17, 0.1 eV 2 < Am 2 < 1 eV 2 . The caveat here is that the CDHS 
limits be applicable in HDM; or, equivalently, that some other experiment (e.g. Super-K) 
limit disappearance in HDM, for the mass splittings of interest for LSND, at least as 
much as CDHS for a given Am 2 . Currently the material in the neutrino path of CDHS 
is unknown. If the experiment had been partially in air, such that Am 2 ^,f S /r > 0.3, 
where r is the fraction of the neutrino pathlength in HDM, and no other experiment places 
as tight of limits on disappearance as the standard CDHS limits, we would have no 
prediction for a signal from such a reactor experiment. 

4.2 2+2 models 

Analysis of the SNO, KamLAND, and Super-K data strongly disfavors 2+2 oscillation 
schemes in the standard interaction picture [35, 37-40] because 2+2 models with LSND 
require large mixing with sterile states, which are disfavored by the SNO, KamLAND and 
Super-K atmospheric data. 

To determine the fit of 2+2 models to the data in HDM, a new analysis would need 
to be carried out. The existing analyses distinguish between active-active versus active- 
sterile oscillations through the MSW effect, which couples only to active neutrinos. These 
analyses, however, are only valid if the MSW effect is the dominant matter effect, and hence 
do not constrain active-sterile mixing in the case that the new non-standard interactions 
considered here dominate. A new analysis with the new interactions, however, would be 
highly model dependent, as the form of the potential due to the new interactions must 
be known. As we have been carrying out only a model independent analysis for such new 
interactions, we will not perform such an analysis here. They could, however, easily be 
carried out once the form of the potential is known. 

4.3 Three neutrino models 

Traditionally, the LSND result has been interpreted to necessitate a fourth, sterile, neu- 
trino because consistency with the Bugey experiment requires /\m\ SND > 0.3 eV, but 
Am 2 tm ~ 3 x 10~ 3 eV, and Am 2 solar ~ 7 x 10 -5 . However, as table [j] indicates, the limits 
from Super-K in HDM are modified so that Am 2 tm < 6 x 10~ 3 eV and the limits from 
LSND so that lS.m 2 LSND > 0.1 eV. This Super-K upward muon data set, however, though 
containing mostly HDM pathlengths, also contains a horizon bin, whose pathlength is 
mostly through air. Without a reanalysis of the data where bins containing pathlengths 
mostly through air are removed, we cannot be sure that the allowed region is an accurate 
reflection of the parameters in HDM. Furthermore, the allowed LSND and Super-K regions 
are approximately an order of magnitude apart. Given such sources for systematic error in 
the analysis and the proximity of the allowed regions, it is interesting to consider models in 
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which both the LSND and Atmospheric results in rock are accommodated with the same 
mass splitting. 

In [9], a model was considered in which the atmospheric and LSND mass splittings were 
accommodated by a splitting of 0(0.1 eV 2 ) in rock, while in air (where the LSND result 
is not applicable) the splitting is O(10 -3 eV 2 ), and the KamLAND result by a splitting of 
0(1(T 4 eV 2 ). In this model, the upper state is described by 

k) = ^(W + K» + e|i/e)) > (4.4) 

to describe the large mixing atmospheric results, with a small admixture of v e for the LSND 
result, and the lower two states by 

and 111 

We discuss here on how such a model would be directly testable. First, in the very 
near future K2K would be able to rule out such a model if it detected high statistics 
disappearance with a mass splitting in the range Am 2 < 0.1 eV 2 . 

The MINOS and OPERA experiments will provide another means to measure this same 
mass splitting in rock in the near future. Additionally, MINOS will have the capability 
of measuring — v e oscillation down to mixings of order, sin 26> Me = 4 x 10~ 2 which arc 
interesting for LSND mass splittings Am 2 < 0.4 eV 2 . 

Like the 3+1 model, this three neutrino model also predicts that a reactor experiment 
with similar sensitivity as Bugey, but with pathlengths instead in HDM, would detect v e 
disappearance in the same parameter range of interest. 



5. Summary 

We have examined the effects on neutrino oscillations of a very light scalar mediating a 
force between neutrinos and ordinary matter, as proposed in [5,6]. As these new matter 
effects could be much larger than the MSW effect, we have considered how these new 
interactions would significantly modify the constraints on the neutrino parameters in any 
High Density Medium (HDM). I summarize in table |2| for each parameter the allowed 
range and experimental signals with standard interactions versus these new in-medium 
non-standard interactions. In particular, indications for these new matter effects from 
already planned and running experiments include: 

• sin 26»g DM with a value lar ger in HDM than the limit given by the CHOOZ experi- 
ment; 

• Am 2 2 , sin 2 2^12 detected by KamLAND in a range not allowed by the SNO results; 

• Am 2 3 , sin 2 2#23 detected by K2K, OPERA, and MINOS with parameters somewhat 
modified from the Super-K result; 
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Model 


Parameter 


SI range 


Experiment 


HDM range 


Experiment 


3 neutrino 


• 2 r\n 

sm 2^i3 





0.1 


1 'I IY \ 1 \ V 





0.5 


buper-K, KamLAIND 
K2K,MINOS,Reactor Expt. 


sm 2^i2 


0.76 


0.88 


blNO 


0.4 


1 


KamLAJND 


sm 2^23 


0.92 


1 


buper-K 


0.8 


1 


buper-K, KzK 
IViilN UD,Ur HjHA 


Am?r, 


6.5 x 10~ 5 


8.2 x 10~ 5 


SNO 


10~ 3 


10~ 5 


KamLAND 


Am| 3 


1.9 x 10~ 3 


3.0 x 10~ 3 


Super-K 


1 x 10" 3 


6 x 10" 3 


Super-K,K2K 
MINOS,OPERA 


3+1 neutrino 


sin 2^ e 


disfavored 


CDHS,Bugey 


10~ 3 


10" 1 


MiniBooNE, Reactor expt. 


2+2 neutrino 


sin 2^ s 


disfavored 


Super-K 





1 


MINOS,OPERA 


sin 29 es 


disfavored 


SNO 





1 


no experiment 



Table 2: Parameters of interest for each type of model. Standard interactions (SI) upper and lower 
limits (90% C.L.) and the experiment which gives that limit (columns three and four). Present 
HDM limits and experiments which arc of interest for measuring these parameters in the near 
future (columns five and six). 



• Detection of v e disappearance in HDM in the region excluded by Bugey, should 
MiniBooNE confirm the LSND signal. 

Any of these signals would be a smoking gun for the new matter effects considered here. 

These upcoming and currently running experiments promise to determine the param- 
eters quite accurately in rock. We have, however, no signal for the parameters in air, as 
well as no information on the form of the density dependence. Without a theory for the 
density dependence of such a force, we are unable to constrain the parameters from solar 
neutrinos, and the constraints from short baseline and atmospheric neutrino experiments 
could be significantly different if the neutrino parameters are sensitively dependent to small 
variations in density. 

Experiments are needed to determine if and how sensitively the neutrino parameters 
vary from such a new force. In the past it has been assumed that new matter effects were 
arising from physics above the standard model energy scale, and for typical earth densities 
were unimportant. It is suggested here that a new attitude be taken in which large non- 
standard in-medium dependences are searched for in terrestrial neutrino experiments. We 
may find that terrestrial neutrino oscillations provide the tool to directly explore a wide 
range of physics beyond the standard model. 
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